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Abstract 

Limits from neutron beta decay on parameters describing physics beyond the Standard 
Model are presented. New Physics is described by the most general Lorentz invariant 
effective Hamiltonian involving vector, scalar and tensor operators and Standard Model 
fields only. Two-parameter fits to the decay parameters measured in free neutron beta 
decay have been done, in some cases indicating rather big dependence of the results on 
9a/ 9v ratio of nucleon form factors at zero four- momentum transfer. 

1 Introduction 

For many years nuclear /3-decays have been exploited as laboratories for testing the Standard 
Model (SM) in the domain of low energies. Along with developments of intense sources of 
cold and ultracold neutrons and improvements of experimental techniques, the precision of 
measurements in the simplest of such systems: the /3-decay of a free neutron, is constantly 
increasing. It opens the way to study the limits on physics beyond SM set solely by the 
parameters of the neutron /3-decay. 

We assume that at the quark-lepton level /3-decay is described by the general 4-point 
Hamiltonian [1] 



k,l = L,R ^ 

+ A kl eP k v {k) uPid 



+ a kk e^=PkV (k) u-^P k d | + H.c, (1) 

where u, d are quark fields, e stands for electron field and Pl = \ (1 — 75), Pr = \ (1 + 75), 
&p,v = § [lfj,,7u], where our metric and gamma matrices are the same as e.g. in [2]. We work 
in the basis in which mass matrix of charged leptons is diagonal and the left (L) and right (R) 
neutrino fields are given by 



(L) 



V 



J2UeiP L Vi, (2) 

i 

E V - P «^' ( 3 ) 



where Vi is the i-th neutrino field with a certain mass, U and V are respectively the Pontecorvo- 
Maki-Nakagawa-Sakata matrix and similar mixing matrix for right-handed neutrinos. SM is 
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restored when aki = Aki = a^k = for k, I = L, R except an = V^Gf/V^, where Gf is the 
usual Fermi constant and V u d is the element of quark Cabibbo-Kobayashi-Maskawa mixing 
matrix. 

When calculating the amplitudes for neutron beta decay at small four-momentum transfer 
q 2 « we have used the relations [1] 

gvUp^u n = (p\u^d\n), (4) 

&4« p 7m75«7i = (p|«7m75^|«} 5 (5) 
9sUpU n = (p\ud\n), (6) 
(p\ua^d\n), (7) 

with (p\, Up and \n), u n being proton and neutron states and bispinors, respectively. From 
conserved vector current hypothesis one gets gy = 1. In the quark model with spherically 
symmetric wave functions of quarks the following relations have been derived [3]: gs = — \ + 
j^gA and gr = § (j% + j^gA^ • Substituting the SM value for g^ — 1.27 into the above relations 
leads to: ~ 0.64 and gx — 1.47. However, in our derivations and fits we treat gs and gx as 
free parameters (independent of 5^4). 



2 Decay Parameters 

From Eq. (1) the five- fold differential decay width for polarized neutron without measurement 
of final electron and proton polarization is given by (in analogy to [4]) 

p e E e E u < 1 + a - + b— 



+ A n 



A^ + B^ + D PeXp 



E e E v E e E u 



(8) 



where A n is the neutron polarization vector, m e , p e = \p e \, E e are, respectively, the mass, 
momentum and total energy of electron, I?o is the maximum value of E e , \p u \ = E v = Eq — E e 
is the antineutrino energy 1 . The Cl e , Q u denotes the solid angles of electron and antineutrino 
emission. We have worked at tree-level (except: calculation of (E^ 1 ) — see below) and with 
approximations such that terms proportional to n p 7su n are not present in (8). Furthermore, 
we will consider only cases when: gy, gA, gs, gr, as well as o^/, A^i, a^k for k, I = L, R are real 
— then D = and time reversal symmetry is preserved, that is well motivated experimentally 
(PDG average [5]: D = (-4± 6) x 10" 4 ). 

We express the decay parameters a, b, A, B, where B has the form of B = Bq + b v m e /E e , 
in terms of the ratio gA/gv and the following parameters (see also [1]) for k, I = L, R 

t t a ki c A k i gs rp ocki gr , n s. 

Vkl — K k , &H — K k , J-kl — K k , (yj 

(ill (ill gv a LL gv 

where 

/t'.iv -i 2 \ 1/2 

n = 1 - w = (iw) ■ (10) 

with summation J2i running only over kinematically allowed antineutrino states. In SM and for 
some cases of physics beyond SM 6 = and b v = 0. As a result of the applied approximations 
formulas for a, b, A, Bq, b u depend in general case only on two combinations 

sl = S LL + S LR , (11) 
sr = S RR + S RL . (12) 



The effect of nonzero neutrino masses enters only trough presence of mixing matrices U and V. 
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Next, following the approach applied in [6, 7, 8], we define 



a{{W- 1 )) 
A({W~ 1 )) 

B{(W- 1 )) 



l + b{W- 1 ) ' 
A 

l + b(W~ l ) 



(13) 
(14) 
(15) 



where (W 1 ) = m e (E e }, and apply these quantities in the fits to the experimental data. The 
X 2 , which we will minimize with the fit procedure, is of the form 



X 



= E 

i 

+ E 

j 

+ E 



a l -a{{W- 1 ). i ) 
5ai 

6A 3 

B k - B{{W~ l ) k ) 



(16) 



where the selected data are presented in Table 1: at, Aj, Bk and <5aj, SAj, 5B k denote the 
central value and the error of the respective decay parameter in a certain experiment. We 
calculate the particular value of (W~ 1 )i = m e (E~ }i from 



dE e 



dT 

dE, 



-El 



jpmax jp 

dE e §r, 

grain CWj e 



(17) 



where EY im and Ef 1 ^ in general are different for different experiments. At this stage of 
calculation Fermi function F(E e ) (that is a leading order QED correction) [9, 10, 11] has been 
incorporated and SM was assumed 



dF 

dEe. 



(9 2 v + 3ffi) ^^ Pe E e (E - E e fF{E e 



F(E e 



2-na E e /p e 
-2iraE e /p e ' 



(18) 
(19) 
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PAR. VALUE ERROR (W 1 ) PAPER ID (PDG) 



a 


-0.1054 


0.0055 


0.655 


BYRNE 


02 


[12] 




-0.1017 


0.0051 


0.655 


STRATOWA 


78 


[13] 

L J 


A 


-0.11966 


0.00166 


0.557 


LIU 


10 


[14] 




-0.1189 


0.0007 


0.534 


ABELE 


02 


[15] 




-0.1160 


0.0015 


0.582 


LIAUD 


97 


[16] 




-0.1135 


0.0014 


0.558 


YEROZOLIMSKY 


97 


[17] 




-0.1146 


0.0019 


0.581 


BOPP 


86 


[18] 


B 


0.980 


0.005 


0.599 


SCHUMANN 


07 


[19] 




0.967 


0.012 


0.600 


KREUZ 


05 


[20] 




0.9801 


0.0046 


0.594 


SEREBROV 


98 


[21] 




0.9894 


0.0083 


0.554 


KUZNETSOV 


95 


[22] 



Table 1: We have followed the PDG [5] data selection but took only the most precise data 
(the error in measurements of a is less than 6% of central value, for A and B — it is less 
than 2%). When experiment report statistic and systematic error separately we add these two 
errors in quadrature. In the case of asymmetric errors we have taken the larger of the reported 
errors. Most of presented values of the {W~ l ) have been taken from [7]. We have used all 
11 "data points" in the table above in every fit presented in this paper. Because of unsolved 
experimental ambiguity of neutron lifetime measurements (see [5]) we have not included this 
quantity in our analyzes. 



3 Results 

In SM the formulas derived for decay parameters depend on A = <7a/sv alone and simplify to 
1 - A 2 2A(1 - A) _ 2A(A + 1) 

A ~ OX2 , 1 > B ~ OX2 i i ■ ( 20 ) 



3A 2 + 1 ' 3A 2 + 1 ' 3A 2 + 1 

In this case, the one-parameter fit is performed, which results in x m i n = 25.42 with 

± 0.0014 (68.27% C.L.) (21) 
A = 1.2703 ± 0.0023 (90% C.L. ) (22) 
± 0.0028 (95.45% C.L.) (23) 

that is in a good agreement with the PDG average [5]: A = 1.2701 ± 0.0025 (error scaled by 
PDG by 1.9). 

The above results (21-23) apply also when: Vki = Ski = Tfcfc = for k, I = L, R except 
Vlr (and Vll — 1 by definition — see Eqs. (9) and (10)). In this case the formulas (20) hold 
for modified A 

A = Mi^k. (24 ) 
9V 1 + V LR 

In the next step one of the parameters: Vr^, Sk, T^k for k = L, R is nonzero and fitted 
together with the ratio qaIqv ■ Among these cases only when the nonzero parameter is sl or 
Tll we have 6/0 and b u ^ 0. The results of such two-parameter fits are presented in the 
Fie. 1. 
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Fi gure 1: Results of the two-parameter fits (the cross denote the position of Xmm ^ 

each 

case). The line — marked areas correspond to the A = qaI 9v intervals: (21) — the narrow 
one and (23) — the wider one. Note that | • | is the absolute value — not the module of a 
complex number, as all parameters are real. 
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In conclusion, Standard Model describes the neutron beta decay very well. The fits are 
minimally better if New Physics is included, especially if tensor terms are present. In some 
cases there is rather big dependence of the results on qaI 9v ratio. 

This work has been supported by the Polish Ministry of Science and Higher Education 
under grant No. N N202 064936. 
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